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I.  INTRODUCTION 


The  effects  of  liquid  payload  motion  on  the  flight  behavior  of  a  pro¬ 
jectile  have  been  of  concern  to  the  Army  for  many  years.  For  the  past  25 
years,  there  has  been  a  systematic  effort  within  the  Army  to  solve  liquid 
payload  problems  and  the  effort  has  included  both  experimental  and  theoretical 
approaches.  Much  of  the  earlier  work  has  been  published  in  Reference  1  which 
is  an  Engineering  Design  Handbook  (1969)  for  liquid  filled  projectiles. 
Reference  2  contains  an  excellent  bibliography  of  much  of  the  significant 
research  including  more  recent  contributions.  Although  many  facets  of  the 
problem  have  been  addressed,  most  of  the  past  effort  has  been  directed  toward 
the  case  of  a  single  component  liquid  payload.  Scott^  considered  a  two- 
component  inviscid  system,  derived  relations  for  the  liquid  eigenfrequencies, 
and  obtained  fair  agreement  with  experiment.  MurphyS  provided  a  formulation 
of  the  effect  of  liquid  motion  on  projectile  stability  by  computing  liquid 
moments  from  the  internal  cavity  wall  pressures  and  wall  shear  stresses.  This 
work  was  applicable  to  a  fully  spun-up  liquid  (solid  body  rotation)  in  a 
completely-  or  partially-filled  cylinder  with  a  solid  or  air  central  core. 
More  recently,  Murphy^  considered  two  immiscible,  viscous  liquids  in  a 
rotating  cylinder.  The  effects  of  viscosity  are  approximated  by  boundary  layer 
treatments  based  upon  the  work  of  Wedemeyer.^  Viscous  effects  are  limited  to 
the  boundary  layer  on  the  cylinder  walls  and  to  the  viscous  shear  layer  at  the 
interface  of  the  two  liquids. 

The  recent  theoretical  advances  for  a  two-component  liquid  and  the  prac¬ 
tical  replications  were  motivating  factors  for  this  experimental  Investiga¬ 
tion.  Therefore,  the  objective  of  this  experimental  investigation  was  to 
provide  an  initial  data  base  for  a  two-component  liquid-filled  cylinder  and  to 
provide  experimental  data  which  could  be  compared  to  theory. 


Engineering  Design  Hmdcook,  Liquid-Filled  Projectile  Design^  M'C  Parphlet 
No.  706-265,  U.S.  Amy  Materiel  Cormand,  Washington,  D.C.,  April  1969. 

(AD  853729) 


Charles  H.  Murphy,  "Angular  Motion  of  a  Spinning  Projectile  with  a  Viscous 
Liquid  Payload,"  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground, 


Maryland,  AF.BP.L-MF.-03294 ,  August  2982.  (AD  A2 25332)  (See  also  Journal  of 
Guidance,  Control,  and  Ditnanics,  Vol.  6,  July-August  2983,  pp.  280-286.) 


3.  W.  E.  Scott,  "The  Inertial  Wave  Speatrun  in  a  Cylindrically  Confined,  In¬ 
viscid,  Incarrpressihle,  Two  Component  Liquid,"  Ballistic  Research  Labora¬ 
tory,  Aberdeen  Proving  Ground,  Maryland,  BRL  Report  1609,  September  1972. 
(AD  752439)  (See  also  Physics  of  Fluids,  Vol.  16,  No.  1,  January  1973, 
pp.  9-12. ) 


4.  Charles  H.  Murphy,  "Side  Moment  Exerted  by  a  Two-Component  Liquid  Payload 
on  a  Spinning  Projectile,"  AIAA  Paper  84-2115,  August  1984.  (See  also 
AEBF.L-TE-02624,  December  1984,  ADA149845.) 


5.  E.  H.  Wedemeyer,  "Viscous  Corrections  to  Stewartson's  Stability  Criterion," 
Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland,  BRL  Reporz 
1325,  June  1966.  (AD  489687) 


7 


11.  EXPERIMENT 


A.  The  Laboratory  Gyroscope. 

The  liquid-filled  cylinder  at  full  spin  is  shown  in  Figure  1.  The  cylin¬ 
der  was  filled  to  90%,  leaving  an  air  core  of  approximately  31%  of  the 
cylinder  radius.  The  liquids  have  different  densities  which  cause  the  heavier 
liquid  to  centrifuge  toward  the  outer  region  of  the  cylinder  and  force  the 
lighter  liquid  inward.  The  two  liquids,  oil  and  water,  are  immiscible  and  an 
interface  is  formed  at  the  radius  bj.  The  cylindrical  container  is  located 

within  the  rotor  of  the  freely  gimballed  gyroscope  shown  in  the  photograph  of 
Figure  2.  The  rotor  is  fairly  massive  and  weighs  approximately  30  times  as 
much  as  liquid.  Rings  of  varying  mass  can  be  added  to  the  rotor  to  change  the 
spin  moment  of  inertia  and,  hence,  change  the  natural  frequency  of  the  system. 
The  rotor  is  mounted  inside  a  non- rotating  cage  with  ball  bearings  at  each 
end.  A  DC,  variable  speed  motor,  attached  to  the  bottom  of  the  cage  drives 
the  rotor  at  the  desired  spin  rate.  A  threaded  spindle  is  situated  on  the  top 
of  the  cage  so  that  non-spinning  weights  can  be  positioned  to  vary  the  trans¬ 
verse  moment  of  inertia,  1^,  which  also  changes  the  natural  frequency  of  the 

system.  The  gimbal  pivots  consist  of  crossed-leaf  springs  (flexural  pivots), 
which  are  instrumented  with  strain  gages.  The  strain  gage  output  is  linear 
with  angle  for  small  yaw  angles.  The  pivots  have  very  small  spring  constants 
and,  therefore,  little  Influence  on  the  gyroscope  motion. 

B.  Data  Acquisition. 

The  axis  of  the  rotor,  at  zero  yaw,  is  in  the  vertical  position.  The 
inner  gimbal,  which  Includes  the  rotor,  is  held  in  the  vertical  position  while 
the  rotor  is  spinning  at  the  desired  rate.  After  sufficient  time  has  elapsed 
for  the  liquid  to  achieve  solid  body  rotation,  the  inner  gimbal  is  released 

with  care  so  as  not  to  impart  a  disturbance.  Any  significant  unstable  motions 
are  self-starting  and  do  not  require  an  initial  disturbance.  As  the  coning 
(yaw)  motion  grows  in  amplitude,  the  yaw  histories  are  recorded  on  a  digital 
oscilloscope.  The  strain  gage  output  from  one  flexural  pivot  is  shown  in 

Figure  3.  If  the  data  are  of  satisfactory  quality,  they  are  stored  on  a 

floppy  disc.  Weights  on  the  spindle  are  then  moved  so  that  another  yaw 

frequency  can  be  obtained  and  the  above  procedure  is  repeated  until  enough 
data  is  acquired  to  define  yaw  growth  rates  over  the  desired  range  of  coning 
frequencies. 

C.  Data  Processing. 

Data  are  transferred  from  the  local  oscilloscope  mediian  to  a  VAX  11-780 
minicomputer  for  processing.  The  data  reduction  program  provides  the  capa¬ 
bility  of  digital  filtering  which  is  usually  necessary;  also,  windowing  is 
sometimes  required  to  discard  obviously  bad  data.  The  program  selects  an 
initial  yaw  amplitude,  Kjq  to  be  used  in  scaling  all  successive  amplitudes, 

Kj.  The  data  are  plotted  as  ln(Kj/Kjg)  versus  time  increment,  and  fitted  by  a 

linear  least  squares  procedure  over  the  desired  Interval.  Figure  4  shows  such 
a  plot  and  the  Interval  of  linearity  Is  seen  to  be  approximately  from  3  to  9 
seconds.  A  linear  variation  of  ln(Kj/Kjg)  with  time  Indicates  a  constant 
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growth  rate.  Usually  the  yaw  ariplitudes  are  small  (in  the  range  of  0.2  to  1.5 
degrees).  Typically,  the  response  of  the  gyroscope  is  carefully  examined  for 

the  maximum  yaw  growth  rate  (tt  )  and  the  coning  frequency  at  which  maximum 

‘max 

yaw  growth  rate  occurs  (td  )•  A  summary  of  the  test  conditions  is  provided 

^ax 

in  Table  1.  Approximately  20  to  24  yaw  histories  were  obtained  for  each  test 
condition. 


TABLE  1.  TEST  CONDITIONS. 


c/a  =  3.126 

a  =  41.21  mm 

I^  =  86.82  g.m" 

f  =  0.9  typical 

Re  =  534000 

value 

^  =  50  rev/sec 

p  =  1.0  g/cc, 

p'  =  .812  g/cc, 
2 

p  =  .960  g/cc, 
2. 

(1.0 

(1.0 

(100 

cs  water) 
cs  oil ) 
cs  oil ) 

Oi 1 /Water 

j  oi  1 

f 

b  /a 

2 

b  /a 

1 

100/0 

1.0  cs 

.895 

.324 

1.000 

85/15 

1.0  cs 

.886 

.337 

.935 

75/25 

1.0  cs 

.886 

.338 

.883 

30/70 

1.0  cs 

.901 

.315 

.607 

15/85 

1.0  cs 

.898 

.319 

.488 

0/100 

1.0  cs 

.901 

.315 

.315 

15/85 

100  cs 

.901 

.315 

.484 

III.  ACCURACY 

Table  2  provides  estimates  for  some  of  the  errors.  The  aspect  ratio 
(c/a)  of  the  cylinder  was  determined  from  physical  measurements  and  is 

believed  to  be  accurate  to  within  ±.05%.  This  small  uncertainty  In  c/a  is 
equivalent  to  almost  1%  error  in  (t  )  and  slightly  more  than  1%  error  in 

^  max 

yaw  growth  rate.  The  fill  ratio,  f,  of  the  cylinder  was  determined  from 

height  measurements  using  a  cathetometer  equipped  with  a  sighting  scope  and  a 

precision  vernier  scale.  The  cathetometer  has  a  resolution  of  .Olmm,  but 
repeated  measurements  indicated  an  uncertainty  on  the  order  of  0.2mm  which 
corresponds  to  less  than  0.1%  in  fill  ratio.  The  yaw  frequency  and  spin  rate 
which  determine  the  dimensionless  yaw  growth  rate  are  believed  to  be  accurate 
and  should  not  be  a  significant  factor  in  evaluating  comparisons  with  theory. 
Yaw  damping  tares,  due  to  the  flexural  pivots,  were  obtained  for  the  empty 
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cylinder  but  were  not  applied  to  the  experimental  data.  !f  the  tare  values 
had  been  applied,  peak  values  of  predicted  growth  rate  would  have  decreased  by 
approximately  1%. 

A  small  offset  of  the  rotor  center  of  gravity  from  the  gimbal  pivot  axis 
produces  a  gravity  moment;  the  restoring  force  of  the  flexural  pivots  also 
produces  a  moment.  These  moments  are  equivalent  to  an  external  moment  which 
influences  the  gyroscopic  stability  of  the  system.  From  some  of  the  tare 
runs,  a  slow  mode  yaw  frequency  could  be  detected  in  addition  to  the  fast  mode 
frequency.  Some  estimates  of  the  gyroscopic  stability  factor  determined  from 
the  two  frequencies  indicated  a  gyroscopic  stability  factor  (s-)  of  almost 
-10.  ^ 


TABLE  2.  ESTIMATED  ERRORS. 


Independent 

Variable 

Error 

Resultant  Error 
in  T,^ 

Resultant  Error 
in  TJ 

c/a 

±.05% 

±.9% 

±1.2% 

f 

±.08% 

±.2% 

±.5% 

• 

±.02% 

±.02% 

coning  rate 

±.03% 

±.03% 

tare  damping 

-.00002 

-1% 

-1/Sg  damping 

-.00006 

-3% 

repeatjbility 

* 

♦Occasional  errors  in  tj  larger  than  the  above  values 


A  value  of  Sg  «  -10  decreases  the  growth  rate  predicted  by  the  theory"*  by 

approximately  3%  but  does  not  noticeably  change  the  yaw  frequency.  The 
computational  and  experimental  data  were  not  corrected  for  the  gyroscopic 
stability  factor.  Computationally,  the  gyroscopic  stability  factor  is  assumed 
to  be  infinite. 

Occasionally  some  values  of  yaw  growth  rate  fall  below  the  well-defined 
trend  of  the  experimental  values.  For  example.  Figure  8  shows  two  obvious 
points  (circles)  which  are  5-10%  below  their  expected  level.  The  cause  of 
this  type  of  error  has  not  been  identified  but  the  number  of  occasions  is;  not 
sufficient  to  alter  the  trend  of  the  data.  None  of  the  identifiable  errors 
are  of  sufficient  magnitude  to  explain  the  observed  frequency  bias  between  the 
experiment  and  theory. 
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IV.  RESULTS  AND  DISCUSSION 


All  of  the  computational  results  shown  in  the  graphical  presentations 
were  obtained  using  Murphy's  theory. The  two-component  liquid  for  most  of 
the  tests  consisted  of  Dow  Corning  1.0  cs  oil  and  water  (also  1.0  cs).  This 
combination  may  at  first  not  seem  too  interesting  because  the  physical  proper¬ 
ties  of  the  liquids  are  not  greatly  different.  However,  computations  show 
that  substantial  changes  in  the  yawing  rates  and  growth  rates  occur  for 
various  oil  to  water  ratios.  Because  of  the  boundary  layer  approximations  in 
theory,  the  low  viscosity  fluids  (which  have  higher  Reynolds  numbers)  should 
give  conditions  where  the  theory  performs  well.  Also,  the  theory,  at  the  time 
of  the  experiments,  was  more  rigorous  for  the  case  of  equal  kinematic  viscosi¬ 
ties.  (Previously,  an  average  Reynolds  number  was  used  for  the  endwall 
boundary  layer  calculations,  but  recent  methods  properly  match  the  viscous 
regions  of  the  two  liquids  on  the  endwalls  of  the  cylinder.) 

The  experimental  results,  in  tabular  form,  are  presented  in  Table  3  (see 
pages  14-17).  The  yaw  growth  rate  parameter  (e),  and  the  fast  mode-  liquid 
side  moment  parameter  (CL5|g{*)are  included  in  the  tabulation  but  were  not  shown 

in  graphical  form  along  with  theory.  Additional  parameters  are  needed  to 
compute  e  and  C^^f^  from  tj  and,  therefore,  decrease  the  accuracy  of  e  ind 

^LSM*  respect  to  tj.  Using  tj  provides  a  more  accurate  comparison  of 

theory  and  experiment.  The  parameters  e  and  C|_5nn*  are  defined  as  follows: 

e  =  *'LSM*  *  e/(’"La2/Ij^) 

D'Amico^  has  shown  that  the  liquid  side  moment  coefficient  is  equivalent  to 
the  above  equation  with  the  following  assumptions:  (1)  no  external  damping  or 
Magnus  moments  exist;  (2)  the  dimensionless  coning  frequency  is  equal  to 
I^/Iy*,  (3)  e  is  relatively  small.  For  the  present  experiment,  values  of 

maximun  CL5f^*  were  typically  6%  smaller  than  predicted  values  of  the  liquid 

side  moment  coefficient. 

Experimental  liquid  eigenf requencies  for  90%  and  100%  filled  cylinder 
obtained  by  Scott3  are  shown  in  Figure  5,  The  parameter  bj/a  indicates  the 

location  of  the  interface  for  the  two  liquids  (the  heavier  fluid  occupies  the 
outer  region  and  the  lighter  fluid  occupies  the  inner  region).  For  bj/a  * 

1.0,  the  cylinder  contains  ail  light  liquid  and  for  bj/a  »  0  (100%  fill)  or 

bj/a  «  .316  (90%  fill)  the  cylinder  contains  all  heavy  liquid.  The  trend  o^ 

the  experimental  results  agree  well  with  theory  but  the  experimental  eigenfre- 
quencies  are  noticeably  higher.  Murphy2  shows  that  a  small  change  in  aspect 


6.  P.  D'Anioo,  ,  "Instabilities  of  a  Gurosaope  Produced  bp  Rapidly 

Rotating ,  Highly  Viscous  Liquids, "  Ballistic  Research  Laboratory ,  Aberdeen. 
Proving  Ground,  Maryland,  AP.BRL~l!R- 02235,  June  19S5.  (AD  A12C374) 
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ratio,  c/a,  of  the  cylinder  can  dramatically  change  the  ya^n  rates.  For  this 
case,  if  the  computational  aspect  ratio  is  Increased  from  3.127  to  3.150  (or 
0.7%),  the  agreement  becomes  very  good  at  all  six  data  points. 

Results  from  the  present  test  are  shown  in  Figures  6  -  Q.  Figure  6  shows 
results  for  three  liquid  compositions  of  oil/water  =  100/0,  85/15,  75/25.  The 
data  illustrate  the  substantial  change  in  growth  rate  and  yawing  frequency 
(t^)  for  the  different  compositions.  The  experimental  data  are  compared  with 

the  computations  which  are  shown  as  the  solid  curves.  The  amplitude  of  the 
growth  rate  agrees  well  with  the  computation  but  the  predicted  yaw  frequencies 
are  consistently  lower  for  all  cases.  The  aspect  ratio  Mas  Increased  by  a 
small  amount  (0,29%)  and  the  computations  were  repeated.  For  the  100/0  case, 
the  correction  was  not  quite  sufficient  and  the  75/25  case  was  slightly  over 
corrected  but  on  the  average  the  bias  was  completely  removed. 

Figure  7  is  another  set  of  data  for  three  liquid  compositions  of  oil/ 
water  *  30/70,  15/85,  and  0/100.  The  results  are  qualitatively  identical  to 
the  results  of  Figure  6.  The  same  consistent  frequency  bias  is  seen  to  exist, 
but  a  0.29%  increase  in  c/a  slightly  overcorrects  the  data  for  all  three 
cases.  The  measured  yaw  frequencies  (t  )  are  about  3%  to  4%  larger  than 

^ax 

predicted  values.  The  identifiable  errors  of  Table  2  are  not  sufficient  to 
account  for  the  discrepancy  but  could  account  for  about  one  third  of  that 
amount.  The  errors  are,  however,  expected  to  be  random  in  nature  and 
therefore  are  not  helpful  in  explaining  the  bias. 

Figure  8  shows  the  results  of  one  experiment  in  which  the  kinematic 
viscosity  of  the  oil  was  substantially  increased.  The  liquid  oil/water  ratio 
was  15/85  with  the  oil  at  a  kinematic  viscosity  of  100  centistokes  (cs)  com¬ 
pared  to  1.0  cs  for  water.  To  illustrate  the  effect  of  the  more  viscous  oil, 
the  experimental  results  for  1.0  cs  oil  are  also  shown  on  the  same  figure.  The 
effect  of  increasing  the  viscosity  is  seen  to  decrease  the  growth  rate  and 
also  to  decrease  t  by  a  substantial  amount.  A  comparison  of  computation 
^ax 

and  experiment  again  shows  a  frequency  bias.  Increasing  the  aspect  ratio  by 
0.29%  brings  the  computed  frequency  up  to  the  experimental  value  but  now  the 
experimental  growth  rate  is  seen  to  be  slightly  low. 

The  values  of  t  were  obtained  at  the  maximum  growth  and  are  shown  in 
Rmax 

Figure  9  as  a  function  of  bj/a  (oil-water  interface).  The  result  is 

qualitatively  similar  to  the  results  of  Scott  (see  Figure  5)  but  the  frequency 
bias  of  the  present  study  is  smaller  by  a  factor  of  2.5  (.29%  vs  .7%)  based  on 
the  c/a  correction.  Figure  9  illustrates  the  precise  agreement  in  the  trend 
of  liquid  eigenfrequencies  with  different  oil/water  compositions.  The  figure 
also  summarizes  the  consistent  frequency  bias  between  experiment  and  computa¬ 
tion. 


The  frequency  bias  illustrated  in  Figures  6  to  9  has  been  previously 
observed  and  is  discussed  in  Reference  2.  The  aspect  ratio  (c/a)  adjustment 
seems  to  correct  the  appropriate  parameters  but  whether  or  not  the  effective 
aspect  ratio  should  be  slightly  different  from  the  measured  value,  has  not  yet 
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been  determined.  A  reasonably  constant  correction  of  0.29"  was  adequate  for 
the  present  experiments  but  other  experiments  have  required  different  correc¬ 
tions. 

Figures  6-8  show  that  an  increase  in  c/a  causes  an  increase  in  the 
predicted  growth  rate  which  could  make  the  experimental  data  appear  low. 
However,  if  the  tare  damping  and  the  gyroscopic  stability  factor,  discussed  in 
the  error  analysis,  were  accounted  for,  the  effect  would  be  to  decrease  the 
predicted  value  and  thus  improve  tho  agreement. 


V.  SUMMARY 

1.  An  initial  set  of  data,  believed  to  he  of  good  quality,  has  been 
acquired  for  gyroscope  experiments  with  two  immiscible  liquids  in  a  right 
ci rcular  cyl  inder. 

2.  The  experimental  data  were  found  to  be  consistent  with  recent 
theoretical  predictions. 

3.  Experimentally  determined  yaw  frequencies  are  consistently  higher 
than  predicted  values  by  three  to  four  percent. 

4.  Further  investigations  are  needed  to  properly  ex.'lain  the  reasons  for 
the  frequency  bias  in  both  single  and  two  component  liquid  systems. 


TABLE  3.  EXPERIMENTAL  DATA 


OIL/WATES  =  100/0,  1.0 

CS  OIL 

OIL/WATER  =  85/15,  1.0 

CS  OIL 

T 

n 

^LSM* 

T 

n 

£ 

^LSM* 

0.0644 

0.000273 

0.00501 

0.1861 

0.0718 

0.000342 

0.00476 

0.1768 

0.0548 

0.000335 

0.00610 

0.2268 

0.0714 

0.000467 

0.00654 

0.2430 

0.0552 

0.000456 

0.00826 

0.3067 

0.0709 

0.000805 

0.01136 

0.4219 

0.0558 

0.000660 

0.01183 

0.4395 

0.0702 

0.000921 

0.01313 

0.4877 

C.0563 

0.000789 

0.01403 

0.5211 

0.0699 

0.001146 

0.01639 

0.6091 

0.0568 

0.003904 

0.01591 

0.5911 

0.0696 

0.001168 

0.01677 

0.6230 

0.0573 

0.000962 

0.01679 

0.6239 

0.0697 

0.001247 

0.01788 

0.6544 

0.0577 

0.000979 

0.01699 

0.6311 

0.0691 

0.001329 

0.01925 

0.7152 

0.0581 

0.000996 

0.01716 

0.6374 

0.0688 

0.001307 

0.01898 

0.7052 

0.0584 

0.000945 

0.01618 

0.6010 

0.0684 

0.001330 

0.01944 

0.7222 

0.0585 

0.000943 

0.01613 

0.5993 

0.0681 

0.001338 

0.01963 

0.7294 

0.0589 

0.000848 

0.01440 

0.5349 

0.0680 

0.001275 

0.01876 

0.6971 

0.0595 

0.000729 

0.01226 

0.4556 

0.0676 

0.001251 

0.01851 

0.6878 

0.0598 

0.000595 

0.00994 

0.3693 

0.0674 

0.001156 

0.01714 

0.6368 

0.0599 

0.000523 

0.00872 

0.3239 

0.0669 

0.001059 

0.01584 

0.5884 

0.0605 

0.000371 

0.00614 

0.2282 

0.0667 

0.000878 

0.01315 

0.4387 

0.0609 

0.000324 

0.00533 

0.1978 

0.0664 

0.000742 

0.01118 

0.4154 

0.0613 

0.000256 

0.00418 

0.1553 

0.0658 

0.000613 

0.00939 

0.348S 

0.0619 

0.000205 

0.00331 

0.1229 

0.0655 

0.000440 

0.00672 

0.2496 

0.0629 

0.000127 

0.00202 

0.0749 

0.0655 

0.000377 

0.00576 

0.2140 

0.0649 

0.000291 

0.00448 

0.1665 

TA?LE  3.  EX-’ERIMENTAL  DATA  (Continued) 


OIL/WATER  =  75 '25,  l.n  i 

CS  OIL 

OIL/WA 

TER  =  30/70,  1.0  1 

CS  OIL 

T 

c 

^LSM* 

T 

e 

^'LSM* 

n.0705 

0.000314 

0.00446 

0.1658 

0.0752 

0.000216 

0.00288 

0.1069 

0.0710 

0.000468 

0.00659 

0.2448 

0.0743 

0.000230 

0.00308 

0.1144 

0.0714 

0.000623 

0.00374 

0.3246 

0.0745 

0.000195 

0.00262 

0.0972 

0.0716 

0.000334 

0.01164 

0.4324 

0.0741 

0.000371 

0.00501 

0.1861 

0.0713 

0.000994 

0.01334 

0.514G 

0.0736 

0.000486 

0.00661 

0.2454 

0.0721 

0.001023 

0.01419 

0.5270 

0.0732 

0.000618 

0.00843 

0.3132 

0.0723 

0.001103 

0.01525 

0.5665 

0.0728 

0.000787 

0.01080 

0.4012 

0.0730 

0.001316 

0.01804 

0.6/02 

0.0726 

0.000960 

0.01323 

0.4916 

0.0742 

0.001389 

0.01871 

0.6952 

0.0724 

0.001040 

0.01436 

0.5335 

0.0731 

0.001403 

0.01918 

0.7126 

0.0721 

0.001102 

0.01528 

0.5673 

0.0745 

0.001434 

0.01924 

0.7149 

0.0720 

0.001229 

0.01707 

0.6341 

0.0745 

0.001510 

0.02027 

0.7529 

0.0718 

0.001323 

0.01842 

0.6844 

0.0739 

0.0U1441 

0.01949 

0.7241 

0.0716 

0.001386 

0.01935 

0.7191 

0.0761 

0.001448 

0.01902 

0.7066 

0.0713 

0.001450 

0.02034 

0.7557 

0.0744 

0.001459 

0.01961 

0.7286 

0.0711 

0.001459 

0.02054 

0.7630 

0.0748 

0.001320 

0.01765 

0.6556 

0.0707 

0.001489 

0.02107 

0.7823 

0.0751 

0.001348 

0.01795 

0.6668 

0.0705 

0.001510 

0.02142 

0.7957 

0.0755 

0.001163 

0.01539 

0.5719 

0.0702 

0.001502 

0.02139 

0.7947 

0.0754 

0.001066 

0.01414 

0.5255 

0.0704 

0.001491 

0.02118 

0.7871 

0.0756 

0.001014 

0.01341 

0.4931 

0.0700 

0.001487 

0.02123 

0.7887 

0.0761 

0.000783 

0.01029 

0.3823 

0.0694 

0.001428 

0.02059 

0.7651 

0.0766 

0.000A36 

0.00569 

0.2113 

0.0690 

0.001371 

0.01936 

0.7378 

0.0686 

0.001204 

0.01756 

0.6523 

0.0680 

0.000914 

0.01343 

0.4990 

0.0668 

0.000331 

0.00570 

0.2120 

TABLE  3.  EXPERIMENTAL  DATA  (Continued) 


OIL/WATER  15/85,  1.0  CS  OIL 


T 

e 

^LSM* 

0.0594 

0.000182 

0.00306 

0.1138 

0.0599 

0.000357 

0.00596 

0.2214 

0.0605 

0.000504 

0.00833 

0.3094 

0.0607 

0.000648 

0.01067 

0.3965 

0.0613 

0.000806 

0.01315 

0.4384 

0.0616 

0.000997 

0.01619 

0.6017 

0.0620 

0.001100 

0.01775 

0.6594 

0.0624 

0.001165 

0.01865 

0.6930 

0.0632 

0.001328 

0.02101 

0.7805 

0.0639 

0.001308 

0.02046 

0.7603 

0.0640 

0.001262 

0.01972 

0.7326 

0.0640 

0.001220 

0.01907 

0.7084 

0.0640 

0.001220 

0.01907 

0.7084 

0.0644 

0.001072 

0.01663 

0.6177 

0.0647 

0.001040 

0.01607 

0.5972 

0.0649 

0.000992 

0.01528 

0.5676 

0.0651 

0.000930 

0.01428 

0.5305 

0.0653 

0.000837 

0.01283 

0.4766 

0.0655 

0.000751 

0.01146 

0.4259 

0.0659 

0.000591 

0.00897 

0.3332 

0.0661 

0.000483 

0.00739 

0.2746 

0.0665 

0.000395 

0.00593 

0.2205 

0.0670 

0.000321 

0.00478 

0.1777 

0.0673 

0.000163 

0.00242 

0.0901 

niL/WATER  =  0/100,  1.0  CS  OIL 


T 

c 

ClSM* 

0.0525 

0.000641 

0.01222 

0.4540 

0.0535 

0.000806 

0.01507 

0.5597 

0.0538 

0.000992 

0.01843 

0.6349 

0.0542 

0.001001 

0.01845 

0.6355 

0.0546 

0.001111 

0.02034 

0.7553 

0.0549 

0.001108 

0.02016 

0.7490 

0.0552 

0.001106 

0.02002 

0.7438 

0.0555 

0.001042 

0.01876 

0.6970 

0.0561 

0.001042 

0.01876 

0.6970 

0.0560 

0.000998 

0.01782 

0.6621 

C.0562 

0.000924 

0.01644 

0.6108 

0.0564 

0.000910 

0.01612 

0.5990 

0.0565 

0.000693 

C. 01226 

0.4557 

0.0568 

0.000710 

0.01251 

0.4648 

0.0569 

0.000592 

0.01041 

0.3867 

0.0570 

0.000445 

0.00781 

0.2902 

0.0573 

n.000486 

0.00849 

0.3155 

0.0575 

0.000335 

0.00582 

0.2162 

0.0579 

0.000342 

0.00591 

0.2194 

0.0581 

0.000320 

0.00551 

0.2046 

0.0584 

0.000210 

0.00360 

0.1336 

0.0585 

0.000200 

0.00341 

0.1268 

0.0594 

0.00009-ir 

0.00158 

0.0535 

TA;^LE  3.  EXPERIMENTAL  DATA  (continued) 


OIL/WATER  =  IS/RS,  IQH  r$  QIL 

T  Tj  e 


0.0517 

0.000320 

0.0524 

0.000419 

0.0531 

0.000543 

0.0535 

0.000599 

0.0540 

0.000674 

0.0545 

0.000726 

0.0549 

0.000724 

0.0552 

0.000700 

0.0554 

0.000758 

0.0557 

0.000758 

0.0560 

0.000691 

0.0563 

0.000742 

0.0568 

0.000672 

0.0571 

0.000604 

0.0578 

0.000524 

0.0585 

0.000446 

0.0592 

0.000370 

0.0599 

0.000274 

0.0607 

0.000181 

0.0613 

0.000184 

0.00620 

0.2303 

0.00800 

0.2971 

0.01033 

0.3837 

0.01119 

0.4159 

0.01249 

0.4640 

0.01332 

0.4948 

0.01319 

0.4901 

0.01268 

0.4710 

0.01367 

0.5079 

0.01361 

0.5053 

0.01234 

0.4585 

0.01318 

0.4898 

0.01183 

0.4397 

0.01057 

0.3925 

0.00907 

0.3370 

0.00762 

0.2832 

0.00625 

0.2320 

0.00457 

0.1700 

0.00298 

0.1108 

0.00300 

0.1114 

Figure  2.  Gyroscope  Photograph. 
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Logarithmic  Plot. 


YAW  GROWTH  RATE 
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Figure  6.  Yaw  Growth  Rates  for  Three  Liquid  Compositions  of 
1.0  cs  on/Water  =  100/0,  85/15,  75/25;  Theory 
(Ref.  4)  and  Experiment. 
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TAU 


Figure  7.  Yaw  Growth  Rates  for  Three  Liquid  Compositions  of 


TAU 


Figure  8.  Yaw  Growth  Rates  for  a  Liquid  Composition  of 
15%  100  cs  Oil  and  85%  1.0  cs  Water;  Theory 
(Ref.  4)  and  Experiment. 
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Figure  9.  Liquid  Eigenfrequencies  for  the  Six  Liquid 
Compositions  of  Figures  6  and  7;  Theory 
(Ref  4)  and  Experiment. 
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LIST  OF  SYMBOLS 


a 

^2 

c 

Clsm* 

f 


•^10 

Re 

C 

V 

"1 

P2 

’^k  ,n 


tT 

^rnax 


^R 


max 


radius  of  a  right  circular  cylindrical  cavity  containing  liquid 

radius  of  the  liquid  interface  for  a  two-component  liquid 

radius  of  the  free  surface  in  a  partially-filled  cylinder 

one-half  of  the  length  of  the  right  circular  cylindrical  cavity 

approximate  liqdd  side  moment  coefficient  for  the  fast  mode 

cavity  fill  ratio;  volume  liquid/volume  cavity 

a^ial  and  transverse  moments  of  inertia  of  the  empty  gyroscope 
rotor 

magnitude  of  the  fast  mode  yaw  arm 
initial  value  of 
Reynolds  number,  a^.t/v 
the  gyroscopic  stability  factor 

kinematic  viscosity 
density  of  the  heavier  liquid 
density  of  the  lighter  liquid 

liq"id  eigenvalue  of  t  for  wave  numbers  k,n;  k  =  l,n  =  3 
for  this  report 

yaw  growth  rate/i),  (the  ordinate  parameter  for  Figures  3, 

4,  and  5 

T|  at  peak  of  bell-shaped  curve 

yawing  rate/J),  (dimensionless  coning  rate) 

tr  at  Tj 

‘max 

rotor  spin  rate 


25 


DISTRIBUTION  LIST 
No.  of 

Copies  Organization 


No.  of 

Copies  Organi zaticn 

12  Administrator 
Defense  Tecnnical 
Information  Center 
ATTN;  OTIC-ODA 
Cameron  Station 
Alexandria,  VA  22304-6145 

1  HODA 

DAMA-ART-M 

Washington,  DC  20310 

1  Commander 

US  Army  Materiel  Command 
ATTN:  AMCDRA-ST 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333-0001 

1  Commander 

Armament  R&D  Center 
US  Army  AMCCOM 
ATTN;  SMCAR-TDC 
Dover,  NO  07801-5001 

1  Commander 

Armament  RSO  Center 
US  Army  AMCCOM 
ATTN:  SMCAR-TSS 
Dover,  NO  07801-5001 

1  Commander 

Armament  R&D  Center 
US  Army  AMCCOM 
ATTN:  SMCAR-LC 
Dover,  NO  07801-5001 

1  Project  Manager 
Armament  R&D  Center 
US  Army  AMCCOM 
ATTN:  AMCPM-CAWS-AM 

Mr.  DellaTerga 
Dover,  NO  07801-5001 

1  0PM  Nuclear 

ATTN:  AMCPM-NUC 

COL.  W.  P.  Farmer 
Dover,  NJ  07801-5001 


4  Commander 

Armament  R&D  Center 
US  Army  AMCCOM 
ATTN*  SMCAR-LCA-F 

Mr.  D.  Mertz 
Mr.  A.  Loeb 
SMCAR-LCA-P 
Mr.  F.  Scerbo 
Mr.  J.  Bera 
Dover,  NJ  07801-5001 

1  Commander 

US  Army  Armament,  Munitions 
and  Chemical  Command 
ATTN:  SMCAR-ESP-L 
Rock  Island,  IL  61299 

1  Oi rector 

Benet  Weapons  laboratory 
Armament  R&D  Centec 
US  Army  AMCCOM 
ATTN;  SMCAR-LCB-TL 
Watervliet,  NY  12189 

1  Commander 

US  Army  Aviation  Research 
and  Development  Command 
ATTN:  AMSAV-E 
4300  Goodfellow  Blvd 
St.  Louis,  MO  63120 

1  Director 

US  Army  Air  Mobility  Research 
and  Development  Laboratory 
Ames  Research  Center 
Moffett  Field,  CA  94035 

1  Commander 

US  Army  Communications  - 
Electronics  Command 
ATTN:  AMSEL-ED 
Fort  Monmouth,  NJ  07703-5301 

1  Commander 

ERADCOM  Technical  Library 
ATTN:  DELSD-L  (Report  Section) 

Fort  Monmouth,  NJ  07703-53C1 


27 


DISTRIBUTION  LIST 


No.  of 


Organization 


No.  of 
Copies 


Organization 


1  AFWL/S!IL 

Kirtland  AFB,  NM  R7117 

1  Conriandant 

US  Army  Infantry  School 
ATTN:  ATSH-CD-CSO-OR 
Fort  Banning,  GA  31905 

1  Commander 

US  Army  Missile  Command 
Research  Development  and 
Engineering  Center 
ATTN:  AMSMI-RD 

Redstone  Arsenal ,  AL  35898-5230 


1  Commandant 

US  Army  Field  Artillery  School 

ATTN:  ATSF-GD 

Fort  Sill,  OK  73503 

1  President 

US  Army  Field  Artillery  Board 
ATTN:  ATZR-BDW 
Fort  Sill,  OK  73503 

1  Commander 

US  Army  Dugway  Proving  Ground 
ATTN:  STEDP-MT 

Mr.  G.  C.  Travers 
Dugway,  UT  84022 


1  Commander 

US  Amy  Missile  Command 
ATTN:  AKSMI-ROK,  Mr.  R.  Deep 

Redstone  Arsenal ,  AL  35898-5230 


1  Commander 

US  Army  Yuma  Proving  Ground 
ATTN:  STEYP-MTW 
Yuma,  AZ  85365 


1  Commander 

US  Army  Tank  Automotive  Command 
ATTN:  AMSTA-TSL 
Warren,  MI  48090 

1  Director 

US  Army  TRADOC  Systems 
Analysis  Activity 
ATTN:  ATAA-SL 
White  Sands  Missile  Range 
NM  88002 


1 


Commander 

US  Army  Development  S  Employment 
Agency 

ATTN:  MODE-TED-SAB 
Fort  Lewis,  WA  98433 


1 


Di rector 

National  Aeronautics  and  Space 
Administration 
Langley  Research  Center 
ATTN:  Tech  Library 
Langley  Station 
Hampton,  VA  23365 


3  Sandia  National  Laboratory 
ATTN:  Mr,  H.  R.  Vaughn 
Dr.  W.  Oberkampf 
Mr.  F.  G.  Blottner 
Albuquerque,  NM  87115 

2  Carco  Electronics 

195  Constitution  Drive 
Menlo  Park,  CA  94025 

1  Commander 

Naval  Surface  Weapons  Center 
ATTN:  Dr.  W.  Yanta 

Aerodynamics  Branch 
K-24,  Building  m-\2 
White  Oak  Laboratory 
Silver  Spring,  MD  20910 

1  Air  Force  Armament  Laboratory 
ATTN:  AFATL/DLODL 
Eglin  AFB,  FL  32542-5000 

1  Aerospace  Corporation 

Aero-Engineering  Subdivision 
ATTN:  Walter  F.  Reddall 
El  Segundo,  CA  9024*= 


28 


No.  of 
Copi es 


DISTRIBUTION  LIST 


Organi zation 


No.  of 
Copies 


Organization 


1  Director 

National  Aeronautics  and  Space 
Administration 

Marshall  Space  Flight  Center 
ATTN:  Dr.  W.  W.  Fowl  is 
Huntsville,  AL  35812 

2  Director 

National  Aeronautics  and  Space 
Administration 
Ames  Research  Center 
ATTN:  Dr.  P.  Kutler 
Dr.  T.  Steger 
Moffett  Field,  CA  94035 

2  Cal  span  Corporation 
ATTN:  G.  Homicz 
M.  Rae 
P.O.  Box  400 
Buffalo,  NY  14225 

1  Raytheon  Company 

Hartwell  Road,  Box  805 
ATTN:  Mr.  V. A. Grosso 
Bedford,  MA  01730 

1  Martin-Marietta  Corporation 
ATTN:  S.H.  Maslen 

1450  S.  Rolling  Road 
Baltimore,  MD  21227 

2  Rockwell  International 

Science  Center 
ATTN:  Dr.  V.  Shankar 

Dr.  S.  Chakravarthy 
P.O.  Box  1085 
Thousand  Oaks,  CA  91360 

1  Arizona  State  University 
Department  of  Mechanical  and 
Energy  Systems  Engineering 
ATTN:  G.P.  Neitzel 
Tempe.  A2  85287 

1  Massachusetts  Institute  of 
Technology 
ATTN:  H.  Greenspan 
77  Massachusetts  Avenue 
Cambridge,  MA  02139 


1  North  Carolina  State  University 
Mechanical  and  Aerospace 
Engineering  Department 
ATTN:  F.F.  DeJarnette 
Raleigh,  NC  27607 

1  Northwestern  University 
Department  of  Engineering 
Science  and  Applied 
Mathematics 
ATTN:  Dr.  S.H.  Davis 
Evanston,  IL  60201 

1  Notre  Dame  University 
Department  of  Aero  Engr 
ATTN:  T.J.  Mueller 
South  Bend,  IN  46556 

1  Rensselaer  Polytechnic 
Institute 

Department  of  Math  Sciences 
Troy,  NY  12181 

1  University  of  Colorado 
Department  of  Astro-Geophysics 
ATTN:  E.R.  Benton 

Boulder,  CO  80302 

2  Univeristy  of  Maryland 
ATTN:  M.  Melnik 

J.D.  Anderson 
College  Park.  MO  20740 

1  University  of  Maryland  - 
Baltimore  County 
Department  of  Mathematics 
ATTN:  Dr.  Y.M.  Lynn 
5401  Wilkens  Avenue 
Baltimore.  MD  21228 

1  University  of  Santa  Clara 
Department  of  Physics 
ATTN :  R,  Greel  ey 
Santa  Clara,  CA  95053 

1  Conmander 

US  Army  Missile  Space 
Intelligence  Center 
ATTN:  AIAMS-YDL 
Redstora  Arsenal,  AL  35898-55C0 


29 


No.  of 
Copies 


DISTRIBUTION  LIST 


No.  of 

Copies  Organization 


Organi zation 


2  University  of  Southern 
Cal i forni a 

Department  of  Aerospace 
Engi neeri ng 
ATTN:  T.  Maxworthy 
P.  Weidman 

Los  Angeles,  CA  90007 

2  University  of  Rochester 
Department  of  Mechanical 
and  Aerospace  Sciences 
ATTN:  R.  Cans 

A.  Clark,  Or. 

Rochester,  NY  14627 

2  Virginia  Polytechnic  Institute 
and  State  University 
Department  of  Aerospace 
Engineering 
ATTN:  Tech  Library 
Dr.  T.  Herbert 
Blacksburg,  VA  24061 

1  University  of  Wyoming 
ATTN:  D.L.  Boyer 
University  Station 
Laramie,  WY  82071 

1  Director 

Lawrence  Livermore  National 
Laboratory 

ATTN:  Mail  Code  L-35 
Mr.  T.  Morgan 
P.O.  Box  808 
Livermore,  CA  94550 

1  University  of  Wisconsin-Madison 
Mathematics  Research  Center 
ATTN:  Dr.  John  Strlckwerda 
610  Walnut  Street 
Madison,  WI  53706 

1  University  of  Virginia 
Department  of  Mechanical 
Ae-ospace  Engineering 
ATTN:  W.  E.  Scott 
Charlottesville,  VA  22904 


10  Central  Intelligence  Agency 
Office  of  Central  Reference 
Dissemination  Branch 
Room  GE-47  HQS 
Washington,  DC  20502 


Aberdeen  Proving  Ground 

Director,  USAMSAA 
ATTN:  AMXSY-D 

AMXSY-MP,  Mr.  H.  Cohen 
AMXSY-R,  Mr.  Robert  Eissner 

Commander,  USATECOM 
ATTN:  AMSTE-TO-F 
AMSTE-CM-F 
Mr.  Gibson  (2  cys) 

PM-SMOKE,  Bldg.  324 
ATTN:  AMCPM-SMK-M 

Cdr,  CRDC,  AMCCOM 
ATTN:  SMCCR-MU 

Mr.  W.  Dee 
Mr.  C.  Hughes 
Mr.  F.  Dagostin 
Mr.  H.  Bach 
Mr.  C.  Jeffers 
Mr.  L.  Shaft 
Mr.  M.  Parker 
ATTN:  SMCCR-RSP-A 

Mr  M.  C.  Miller  (2  cys) 
ATTN:  SMCCR-SPS-IL 
SMCCR-RSP-A 
SMCCR-MU 


30 


V  '"L*^  *.•  O  '^'  •■■’  ^  '^•l  '^«  • 


•  -fmj-  gm  ;•  '-.•  •'.• 


USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 

This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the 
reports  it  publishes.  Your  comments/answers  to  the  items/questions  below  will 
aid  us  in  our  efforts. 

1.  BRL  Report  Number _ Date  of  Report _ 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  reed?  (Comment  on  purpose,  related  project,  or 

other  area  of  interest  for  which  the  report  will  be  used.) _ 


4.  How  specifically,  is  the  report  being  used?  (Information  source,  design 
data,  procedure,  source  of  ideas,  etc.) _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far 
as  man-hours  or  dollars  saved,  operating  costs  avoided  or  efficiencies  achieved, 
etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future 
reports?  (Indicate  changes  to  organization,  technical  content,  format,  etc.) 


i 


I 


Name 


CURRENT 

ADDRESS 


Organization 


Address 


City,  State,  Zip 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the 
New  or  Correct  Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


I 

* 


OLD 

ADDRESS 


Name 


Organization 


Address 


City,  State,  Zip 


(Remove  this  sheet  along  the  perforation,  fold  as  indicated,  staple  or  tape 
closed,  and  mail.) 


I 

i 


I 

« 

! 


■  Director 

U.S.  Army  Ballistic  Research  Laboratory 
ATTN;  SLCBR-DD-T 

I  Aberdeen  Proving  Ground,  MD  21005-5066 

I  _ 

OTFICIAL  BUSINESS 
nSNALTY  rON  MIVATC  use. 


I 

I  Director 

’  U.S.  Army  Ballistic  Research  Laboratory 

i  ATTN:  SLCBR-DD-T 

'  Aberdeen  Proving  Ground,  MD  21005-9989 


BUSINESS  REPLY  MAIL  I 

POST  CLASS  PERMIT  NO  12062  WaSHlWCTOW,0C  | 
postage  will  BC  paid  Bt  OePftRTMENT  OP  THE  ARMY 


I 


I 


I 


